Optical switching of radical pair conformation enhances magnetic sensitivity 
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The yield of chemical reactions involving intermediate radical pairs is influenced by magnetic fields 
well beyond the levels expected from energy considerations. This dependence can be traced back to 
the microscopic dynamics of electron spins and constitutes the basis of the chemical compass. Here 
we propose a new experimental approach based on molecular photoswitches to achieve additional 
control on the chemical reaction and to allow short- time resolution of the spin dynamics. Our 
proposal enables experiments to test some of the standard assumptions of the radical pair model 
and improves the sensitivity of chemical magnetometers by two orders of magnitude. 



Introduction. — Radical pair (RP) reactions have re- 
cently raised attention as one of the best characterized in- 
stances in the field of quantum biology. Due to the strong 
dependence of the products' chemical yield on the exter- 
nal magnetic field, they constitute one of the two main 
hypotheses to explain how certain species of birds use the 
Earth's magnetic field to orient during migration. Recent 
advancements are the identification of a suitable protein, 
which is present in the retina and possesses the proper- 
ties required by the avian magneto-reception [TJ [2] , and 
the experimental demonstration that a RP reaction can 
be sensitive to magnetic fields as weak as the Earth's one 
[3]. The formulation of an additional quantum physics 
perspective boosted the interest towards the study of en- 
tanglement, its decoherence and the implementation of 
quantum control techniques [4}{7]. In this respect, a RP 
reaction constitutes an electron spin system exposed to 
a mesoscopic, non-Markovian nuclear spin bath. 

The theory describing the radical pair chemical reac- 
tions has been developed in the spin chemistry commu- 
nity and has successfully undergone a series of experi- 
mental tests [8, 9 . Still, it would be necessary to inves- 
tigate the model on the short timescale of the electron 
spin dynamics. Here, we propose a new experimental 
setup which allows sub-nanosecond time resolution by 
effectively engineering the reaction kinetics. We suggest 
to connect the radicals via a molecular switch presenting 
two isomeric forms and capable of undergoing conforma- 
tional changes upon the absorption of a photon. In this 
way, one can optically manipulate the distance between 
the radicals and effectively determine the instant of their 
re-encounter. Since the chemical reaction takes place at 
the moment of re-encounter, we achieve control on the 
reaction kinetics. This new possibility adds in a comple- 
mentary way to previous proposals in which the electron 
spin dynamics is influenced by the application of mag- 
netic pulses [4 or oscillating fields 0. 

First, we present the experimental setup and discuss 
its feasibility with state-of-the-art technology. The con- 
sistency of time scales and the efficiency of the procedure 
are analyzed introducing the case study of azobenzene as 
the photo-controlled bridge. Second, we show how the 



additional control can be exploited to improve the per- 
formances of sensors based on the RP mechanism, and 
in particular we apply our scheme to a chemical mag- 
netometer to boost its sensitivity by two orders of mag- 
nitude. It is remarkable that such enhancement is even 
obtained using a conservative, i.e. coarse-grained, time 
resolved scheme. 

Radical Pair Chemical reactions. — The name RP re- 
actions refers to a broad class of chemical reactions 
characterized by the creation of radicals as intermediate 
products [10 . Here, we consider those reactions which 
take place in solution and involve two different chem- 
ical compounds, let us call them A and D for accep- 
tor and donor respectively, free to diffuse. Typically, 
one of the molecules is photoexcited and, if A and D 
are close enough, a fast electron transfer establishes the 
radical pair D* + -A*~. The electron spins are assumed 
to be in the singlet state \S) = ^(|tl) - lit)), while, 
at room temperature, the nuclear spins are in the com- 
pletely mixed state. The radicals diffuse away from each 
other to such distances that the electron-electron inter- 
actions become negligible and the electron spin dynamics 
is determined only by the external magnetic field and the 
interaction of each electron with the surrounding nuclear 
spins on the respective radical [TO] : 

k=l,2 k k,j 

where j e = —g e ^B is the electron gyromagnetic ratio, 
Sk, Ikj are the electron and nuclear spin operators re- 
spectively, B is the external magnetic field, c\\ Zj denote 
the hyperfine coupling tensors, and we fix h = 1. In a 
stochastic way, geminate re-encounter of the two radi- 
cals happens, and backward electron transfer completes 
the chemical reaction for those radicals in a singlet state. 
Triplet radicals recombine to distinct triplet products. 

The number of radical pairs which recombine through 
the singlet reaction channel is directly proportional to 
the time integral of the fraction of radicals in a sin- 
glet state f s (t) weighted by the probability of (diffusive) 
re-encounter: <I> S = J °° p re (t)f s (t)dt , where f s (t) = 
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(S\ Pei(t) \S) is the overlap between the singlet state and 
the electron spin state p e i(t) at time t. Thus, we are 
facing a system in which two distinct, but concurrent 
dynamics are taking place: The quantum evolution of 
the electron spin state, and the classical diffusion of the 
molecules D ,+ and A* - . The latter one stochastically 
determines the duration of the former. In the literature, 
the probability distribution for the re-encounter time has 
phenomenologically been described as a single exponen- 
tial p re (t) = ke~ kt [10j[TT] with re-encounter rate k. It is 
Pre(t) that we aim to control in order to modify <£ s and 
all the related quantities. 

Re-encounter Probability. — A large rate k ~ Ins -1 
indicates that most of the geminate re-encounters hap- 
pen in the first few nanoseconds after the initial electron 
transfer, and so it becomes hard to study the electron 
spin dynamics at later times or to observe the depen- 
dence of the entanglement lifetime on B as predicted in 
[4]. Furthermore, a weak magnetic field effectively has 
no time to produce any appreciable effect and this limits 
the sensitivity of chemical magnetometers or compasses. 
In the context of avian magneto-reception, the RP is 
actually required [TJ [12] to have a small re-encounter 
rate k ~ 1/is -1 . Unfortunately this is generally not 
the case for chemical reactions involving bi-ionic radi- 
cals in solution since the Coulomb attraction leads to 
a fast re-encounter. Thus several experimental efforts 
have been undertaken to modify the distribution of the 
RP re-encounter time, mainly by imposing constraints on 
the geometry of the system as shown in Fig. [I] (a) The 
radicals are linked by a flexible, inextensible chain whose 
length determines the maximum separation that they can 
achieve [13] ; (b) The reaction takes place inside a micelle 
which effectively defines a much smaller container where 
the radicals are allowed to diffuse [14 . The experimen- 
tal studies of RP reactions in these systems confirm a 
substantial modification of the re-encounter time distri- 
bution, but are still interpreted within the phenomeno- 




( a) flexible chain (b) micelle 




(c)photoswitching bridge 

Figure 1: Three approaches to modify the re-encounter prob- 
ability: (a) The radicals are linked by a flexible, inextensible 
chain, (b) the chemical reaction happens inside a micelle, (c) 
the radicals are linked by a rigid, optically controlled bridge. 



logical exponential model with rate k ~ 0.1 — 1 /is -1 es- 
sentially determined by solvent viscosity, chain length, 
hydrophobic tail of the surfactant, and size of the mi- 
celles. These properties cannot in general be changed 
continuously or even independently, and their effect on 
the rate k is not easy to predict quantitatively. To over- 
come these difficulties, we propose to use optical switches, 
in this way combining geometrical constraints with the 
possibility of optically controlling them. 

Photo- controlled radicals. — The distance of the radi- 
cals strongly affects both the possibility of recombination 
and the electron spin dynamics. At small separations, the 
electron transfer is favoured while the electron spin dy- 
namics is suppressed due to the large exchange and dipo- 
lar interaction; conversely, at large distances, the electron 
transfer is inhibited while the dynamics is determined by 
the Zeeman and hyperfine interaction only, since the ex- 
change and dipolar interactions are now negligibly weak. 
We propose to optically control the inter-radical separa- 
tion by chemically attaching A and D to the two end- 
ings of a molecule exhibiting two isomeric forms and an 
isomerization pathway which is photo-activated. Such 
molecules possess the property of reversibly switching 
between a straight and a contracted structure upon the 
absorption of a photon of suitable wavelength. The two 
isomeric forms effectively impose a "closed" , respectively 
"open", configuration for the radicals, see Fig. [I] that we 
are able to change in a reversible way by means of short 
laser pulses. Examples of photoswitchable compounds 
are fulgides, diarylethenes or azobenzene derivatives [15] . 

For this description to hold, it is necessary that the 
photo-isomerization process is effective and causes a sub- 
stantial change of the distance between the two radicals. 
In this case, and via the optical manipulation of the open- 
ing and closing of the linking bridge, it is in principle pos- 
sible to directly control the instant of RP creation, the 
exact duration of the quantum dynamics for the elec- 
tron spins, and the instant of the final recombination. 
The implementation of such a scheme is ideally as fol- 
lows: Initially all the bridges are closed and the radicals 
not yet formed, then the RP are created by excitation 
via a laser pulse (standard technique applied for exam- 
ple in [3 ), and immediately after that, a second laser 
pulse triggers the isomerization and separates the rad- 
icals. The electron spin state evolves according to the 
Hamiltonian in Eq. ([!]) until a third laser pulse closes 
the bridges and leads to the final recombination. It is 
desirable, but not strictly necessary, that all the three 
processes could be selectively triggered by laser pulses of 
different wavelengths. 

Under the condition that the duration of the various 
pulses remains much shorter than the time scale of the 
electron spin evolution, the re-encounter probability can 
be described by a peaked distribution p re (t) oc 5(t — r re ), 
with r re being the time elapsed between the second and 
third pulse and "oc" indicating the proportionality with 
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respect to the actual photo-isomerization yield. Such 
peaked distribution allows to monitor the electron spin 
dynamics at the precise instant r re which can be cho- 
sen without restrictions in a wide interval ranging from 
a few nanoseconds to several tens of microseconds (for 
later times one has to take into account spin-lattice re- 
laxation) . The proposed scheme can be straightforwardly 
adapted to engineer more complicated and structured re- 
encounter probability distributions: A series of pulses 
creates a re-encounter probability exhibiting separated 
peaks, while continuous illumination gives rise to an ex- 
ponential distribution (here not phenomenologically as- 
sumed, but externally imposed) whose rate is directly 
proportional to the light intensity. Combining laser 
pulses, dark intervals, and continuous illumination with 
adjustable light intensity, the ability of engineering a spe- 
cific p re (t) is practically limited by technical restrictions 
(e.g. the laser intensity or the speed of light modulation) 
and by the absorption and isomerization properties of the 
photoswitches (e.g. the selectivity with which an "open- 
ing pulse" does not cause spurious "closing events" , and 
viceversa). Realistically, a fraction of the photoswitches 
may not undergo phot oisomerizat ion, but, for thermally 
stable isomers, this means that the corresponding radi- 
cals simply do not contribute to the chemical reaction. 
The actual isomerization yield determines the number of 
radical pairs effectively involved in the experiment. 

Note, that an alternative scenario is given by pho- 
toswitches that do not change their conformation, but 
rather change their conductance, passing from an insu- 
lating to a conductive state and viceversa [I6j [17] . In the 
case of backward electron transfer through the linker, 
the isomerization of the linker corresponds to switching 
on/off the coupling between the two radicals at its ends. 

Case study: azobenzene. — To explore the actual feasi- 
bility of our proposal, we consider the concrete example 
of azobenzene as the photoswitching molecular bridge. 
Azobenzene is an organic molecule composed of two 
phenyl rings linked by a N=N double bond, and occurs 
in two isomeric forms, trans and cis azobenzene, having 
a stretched, respectively contracted, structure. The syn- 
thesis of azobenzene derivatives having different chemi- 
cal functional groups extending from the phenyl rings is 
routinely performed p~8j [19] . Preliminary checks have to 
verify that the chemistry of the radical pair is not drasti- 
cally altered by the presence of the azobenzene and that 
the constraints on the RP distance are satisfied. 

For a realistic implementation of our proposal, the in- 
volved processes need to satisfy a hierarchy of time scales: 
The electron spin dynamics determines the reference time 
scale to be ~ 1 — 10 ns (for typical values of the external 
magnetic field and hyperfine couplings, i.e. ~ 0.1 — 1 mT) 
with respect to which the spin lattice relaxation and the 
thermally driven isomerization have to be slow, while the 
radical creation and photoisomerization mechanism have 
to be fast. Typically, the spin-lattice relaxation takes 



~ 1 — 10 /is at room temperature, while the donor photo- 
excitation and successive electron transfer happens in a 
few picoseconds. In the case of azobenzene, the isomers 
can be considered completely stable at room tempera- 
ture on the time scale of seconds [20] , while experimental 
studies with azobenzene in solution afforded Si excited 
state life times of trans- azobenzene of ~ 2.5 ps when ex- 
cited without excess vibrational energy [21] and ~ 0.5 ps 
when from higher vibrational states of Si populated after 
internal conversion from S2 [22] . Finally, the time scale 
over which we can modulate the re-encounter probability 
depends on the duration of the optical pulses: The esti- 
mates in [23] suggest that a single pulse of 150 fs FWHM 
at 439 nm and having an energy of 200 nJ leads to 3% 
cis-to-trans photoisomerization yield. A single pulse of 
duration ~ 10 ps and analogous intensity can achieve 
photoisomerization yield of 30 — 40%. Concerning the 
process of RP creation, in reference [3] radicals were gen- 
erated by a 7 ns , 5mJ pulse, but even radical pair gen- 
eration in less than one picosecond has been achieved in 
experiments with femtosecond laser pulses [24] . 

Application: chemical magnetometry. — The control of 
the RP re-encounter probability finds immediate appli- 
cation to increase the performance of chemical devices. 
Here, we show how a simple-to-implement control scheme 
highly enhances the sensitivity of a chemical magnetome- 
ter by two orders of magnitude. The basic idea behind 
the functioning of a chemical magnetometer is that, since 
a change in the magnetic field modifies the amount of 
singlet products, one can reverse the reasoning and mea- 
sure the chemical yield to estimate B. By definition, the 
magnetic sensitivity is high when a small change in the 
magnetic field intensity produces large effects on the sin- 
glet yield: 

K{B) = —j±-L= J Pre (t)g s (B,t)dt, (2) 

with g s (B,t) = d ^ S QB^ being the instantaneous mag- 
netic sensitivity. The functional form of f s (B,t) = 
(S\ pei(t) \S) strongly depends on the specific realization 
of the radical pair, in particular on the number of the 
surrounding nuclear spins. Here, we consider a radical 
pair in which the first electron spin is void of hyper- 
fine interactions, while the second electron spin inter- 
acts isotropically with one spin-1 nucleus. In the context 
of the chemical compass (i.e. when the task is deter- 
mining the magnetic field direction through anisotropic 
hyperfine interactions), an analogous configuration (with 
only one spin- 1/2 nucleus) has been proposed [2 , and 
numerically characterized [7 , as being optimal: Addi- 
tional nuclear spins would perturb the intuitive "refer- 
ence and probe" picture. The Hamiltonian then simpli- 
fies to H = ->y e B(s[ z) + S { 2 z) ) + |7e|a5 2 • /, where a is 
the isotropic hyperfine coupling. 

We evaluate the sensitivity in the regime of long spin- 
lattice relaxation times and weak magnetic fields, i.e. 
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Figure 2: Instantaneous magnetic sensitivity for a — 1 mT 
and B = The peaks of the envelope grow linearly with 
time and the envelope's shape is approximately periodic with 
period r en vei = j^Jb ~ 1-07ms. Inside the envelope, g s (t) 
oscillates with period ~ 3p^ru ~ 23.8ns (see inset). 



for 5 <C a, noting that in practical applications one 
can reproduce such condition by applying a suitable ex- 
ternal field for compensation. For the realistic choice 
a 1.0 mT, the Earth magnetic field B Earth — ot/20 is 
indeed weak. Straightforward, but rather tedious, calcu- 
lations lead to an analytical expression for the instanta- 
neous magnetic sensitivity g s (B 1 t) [25 , whose behavior, 
in the regime is shown in Fig. |2j 

When integrated on timescales longer than the effec- 
tive "period" of the envelope, r enve i = ■ 3 T R , the positive 

ITe 1-° 

and negative contributions cancel each other and there- 
fore drastically reduce the ultimate sensitivity A S (B). A 
possible solution is to engineer the re-encounter probabil- 
ity in such a way that the radical pairs are allowed to re- 
combine only when g s (B,t) > 0. Such desirable p re (t) is 
obtained applying weak light pulses (for the final step of 
"closing" the photoswitchable bridge) interspersed with 
equally long dark intervals, both of duration r ~ r enve i/2. 
Considering the realistic parameters of Fig.[2j one obtains 
r « 0.54 /is corresponding to a laser repetition of around 
1 MHz. The corresponding re-encounter probability is 
a piecewise exponentially decreasing function with decay 
rate k prot directly proportional to the light intensity. The 
phenomenological and controlled re-encounter probabili- 
ties are shown in Fig. [3] together with the corresponding 
integrand involved in the evaluation of the magnetic sen- 
sitivity via Eq. 

The change in the (integrated) magnetic sensitivity is 
very sensitive to the duration of the pulses, so we propose 
to monitor A s (B,r) (as currently done in spin chemistry 
experiments [26, 27 ) while changing the pulse duration 
r. From an irregular landscape, a sharp peak centered 
around a certain r' emerges due to the resonant condi- 
tion t' ~ T enve i/2, see Fig. [4] (left) in which r' = 0.535 /is. 
As expected, at r c± 3r' , a weaker resonance gives rise to 
a second peak. Thus, a first estimate of the field inten- 
sity is given by B est im = 217 V 7 ' anc ^ ^ S P rec ^ se vame is 
obtained by applying pulses of duration r' . Fig. [4] (right) 
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Figure 3: Re-encounter probability (top) and integrand for 
the magnetic sensitivity (bottom) over time. The blue 
curve (left) corresponds to the phenomenological exponen- 
tial model, the red curve (right) to the optically controlled 
case. Note that we have chosen the laser intensity such that 
kprot — k — 1/is -1 to make the comparison more direct. In 
general, while k is determined by the radicals and properties 
of the solution, k pro t is a tunable parameter of the protocol. 



shows that the optical control of the RP re-encounter 
probability enhances the magnetic sensitivity up to two 
orders of magnitude with respect to the phenomenolog- 
ical exponential model over a broad range of rate k. A 
similar enhancement is obtained for several k prot satisfy- 
ing kp ro t < r _1 . Similar results are obtained for a single 
nuclear spin-1/2 and a suitably modified protocol [25] . 

Summary and Outlook. — We have proposed a new ex- 
perimental approach which enables the control of the re- 
encounter probability of radical pairs in solution via light 
pulses. This will allow one to study the electron spin 
dynamics on the short time scale of its evolution and 
thus enlarge the information accessible in spin chemistry 
experiments. The control of the reaction kinetics con- 
stitutes a complementary and autonomous method com- 
pared to previous works proposing direct electron spin 
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Figure 4: Left: Magnetic sensitivity A s vs. pulse duration r. 
B = 0.05 mT and l/k pro t = 2 /is (solid), 1 /is (dashed), 0.5 /is 
(dotted). Right: Magnetic sensitivity A s vs. field intensity 
B for the exponential model with 1/k = 10ns (solid, blue), 
1/k — 100ns (dashed, blue), 1/k — 2 /is (dotted, blue), and 
for our proposal with pulse duration r' = 0.535 /is (solid, red). 
Parameters as in Fig. [2] 
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manipulation [H |6j [7] . Such an additional handle on the 
RP dynamics should contribute to a better understand- 
ing of the radical pair mechanism itself. 

In order to improve the performance of man-made sen- 
sors based on chemical reactions, one naturally pursues 
the approach of artificially synthesizing suitable com- 
pounds like, in our proposal, radical pair molecules con- 
nected to a photoswitchable bridge. As a concrete re- 
sult, we have presented a simple protocol which enhances 
the sensitivity of a chemical magnetometer by up to two 
orders of magnitude. Remarkably, this scheme can be 
implemented using weak laser pulses having a duration 
much longer than the time scale given by the electron 
spin dynamics. 
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